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We  analyze  single-shot  readout  for  superconducting  qubits  via  controlled  catch,  dispersion,  and  release  of 
a  microwave  field.  A  tunable  coupler  is  used  to  decouple  the  microwave  resonator  from  the  transmission  line 
during  the  dispersive  qubit-resonator  interaction,  thus  circumventing  damping  from  the  Purcell  effect.  We  show 
that  if  the  qubit  frequency  tuning  is  sufficiently  adiabatic,  a  fast  high-fidelity  qubit  readout  is  possible  even  in  the 
strongly  nonlinear  dispersive  regime.  Interestingly,  the  Jaynes-Cummings  nonlinearity  leads  to  the  quadrature 
squeezing  of  the  resonator  field  below  the  standard  quantum  limit,  resulting  in  a  significant  decrease  of  the 
measurement  error. 


PACS  numbers:  03.67.Lx,  03.65.Yz,  42.50.Pq,  85.25.Cp 

Introduction. -Fast  high-fidelity  qubit  readout  plays  an  im¬ 
portant  role  in  quantum  information  processing.  For  super¬ 
conducting  qubits  various  nonlinear  processes  have  been  used 
to  realize  a  single-shot  readout  [1-6].  Linear  dispersive  read¬ 
out  in  the  circuit  quantum  electrodynamics  (cQED)  setup 
[7,  8]  became  sufficiently  sensitive  for  the  single-shot  qubit 
measurement  only  recently  [9,  10],  with  development  of  near- 
quantum-limited  superconducting  parametric  amplifiers  [9- 
11].  In  particular,  readout  fidelity  of  94%  for  flux  qubits  [9] 
and  97%  for  transmon  qubits  [10]  has  been  realized  (see  also 
[12]).  With  increasing  coherence  time  of  superconducting 
qubits  into  10-100  ps  range  [13,  14],  fast  high-fidelity  readout 
becomes  practically  important,  for  example,  for  reaching  the 
threshold  of  quantum  error  correction  codes  [15],  for  which 
the  desired  readout  time  is  less  than  100  ns,  with  fidelity  above 
99%. 

A  significant  source  of  error  in  the  currently  available 
cQED  readout  schemes  is  the  Purcell  effect  [16]  —  the  cavity- 
induced  relaxation  of  the  qubit  due  to  the  always-on  coupling 
between  the  resonator  and  the  outgoing  transmission  line.  The 
Purcell  effect  can  be  reduced  by  increasing  the  qubit-resonator 
detuning;  however,  this  reduces  the  dispersive  interaction  and 
increases  measurement  time.  Several  proposals  to  overcome 
the  Purcell  effect  have  been  put  forward,  including  the  use  of 
the  Purcell  filter  [17]  and  the  use  of  a  Purcell-protected  qubit 
[18].  Here  we  propose  and  analyze  a  cQED  scheme  which 
avoids  the  Purcell  effect  altogether  by  decoupling  the  res¬ 
onator  from  the  transmission  line  during  the  dispersive  qubit- 
resonator  interaction. 

Main  idea  and  results-  Similar  to  the  standard  cQED  mea¬ 
surement  [7-10],  in  our  method  (Fig.  1)  the  qubit  state  af¬ 
fects  the  dispersive  shift  of  the  resonator  frequency,  that  in 
turn  changes  the  phase  of  the  microwave  field  in  the  resonator, 
which  is  then  measured  via  homodyne  detection.  However, 
instead  of  measuring  continuously,  we  perform  a  sequence  of 
three  operations:  “catch”,  “disperse”,  and  “release”  of  the  mi¬ 
crowave  field.  During  the  first  two  stages  a  tunable  coupler 
[19,  20]  decouples  the  outgoing  transmission  line  from  the 
resonator.  This  automatically  eliminates  the  problems  associ¬ 


FIG.  1 :  (color  online),  a)  Schematic  of  the  measurement  setup.  The 
radio  frequency  (RF)  source  produces  a  microwave  pulse,  which  pop¬ 
ulates  the  resonator  via  a  small  capacitor  Cm.  The  resonator  photons 
then  interacts  with  a  capacitively  (Cg)  coupled  qubit.  The  interac¬ 
tion  with  the  outgoing  transmission  line  is  controlled  by  a  tunable 
coupler,  which  releases  photons  at  the  end  of  the  procedure.  The  re¬ 
leased  field  is  then  amplified  and  mixed  with  the  local  oscillator  (LO) 
signal  to  be  measured  via  homodyne  detection,  b)  The  RF  pulse  B(t) 
(blue  curve)  and  varying  qubit  frequency  cjq(i)  (red  curve),  with  ap¬ 
proximate  indication  of  the  “catch”,  “disperse”,  and  “release”  stages. 
Dashed  lines  show  the  resonator  frequency  uir  and  initial/final  qubit 
frequency  cjq;  A  =  cur  —  coq  is  the  detuning  at  the  “disperse”  stage. 


ated  with  the  Purcell  effect,  as  coupling  to  the  incoming  mi¬ 
crowave  line  can  be  made  very  small  [20]. 

During  the  “catch”  stage,  the  initially  empty  resonator  is 
driven  by  a  microwave  pulse  and  populated  with  ~10  pho¬ 
tons.  At  this  stage  the  qubit  is  far  detuned  from  the  resonator 
[Fig.  1(b)],  which  makes  the  dispersive  coupling  negligible 
and  allows  the  creation  of  an  almost-perfect  coherent  state  in 
the  resonator.  At  the  next  “disperse”  stage  of  the  measure¬ 
ment,  the  qubit  frequency  is  adiabatically  tuned  closer  to  the 
resonator  frequency  to  produce  a  strong  qubit-resonator  in- 
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teraction  (it  may  even  be  pushed  into  the  nonlinear  regime). 
During  this  interaction,  the  resonator  field  amplitudes  (Aeff) 
associated  with  the  initial  qubit  states  |0)  and  |1)  rapidly  ac¬ 
cumulate  additional  phases  and  separate  in  the  complex  phase 
plane  [see  Fig.  2(a)].  Finally,  at  the  last  “release”  stage  of 
the  measurement,  after  the  qubit  frequency  is  again  detuned 
from  the  resonator,  the  resonator  photons  are  released  into  the 
outgoing  transmission  line.  The  signal  is  subsequently  ampli¬ 
fied  (by  a  phase-sensitive  parametric  amplifier)  and  sent  to  the 
mixer  where  the  homodyne  detection  is  performed. 

With  realistic  parameters  for  superconducting  qubit  tech¬ 
nology,  we  numerically  show  that  the  measurement  of  30-40 
ns  duration  can  be  realized  with  an  error  below  10~3,  neglect¬ 
ing  the  intrinsic  qubit  decoherence.  The  latter  assumption  re¬ 
quires  the  qubit  coherence  time  to  be  over  40  /is,  which  is 
already  possible  experimentally  [14].  It  is  interesting  that  be¬ 
cause  of  the  interaction  nonlinearity  [21,  22],  increasing  the 
microwave  field  beyond  ~  10  photons  only  slightly  reduces  the 
measurement  time.  The  nonlinearity  also  gives  rise  to  about 
~50%  squeezing  of  the  microwave  field,  which  provides  an 
order-of-magnitude  reduction  of  the  measurement  error. 

The  model. -We  consider  a  superconducting  phase  or  trans¬ 
ition  qubit  capacitively  coupled  to  a  microwave  resonator  [Fig. 
1(a)].  For  simplicity  we  start  with  considering  a  two-level 
qubit  (the  third  level  will  be  included  later)  and  describe  the 
system  by  the  Jaynes-Cummings  (JC)  Hamiltonian  [7]  with  a 
microwave  drive  (we  use  h  =  1) 

H  =  wq(f)cr+(T_  +  u>Ta'a  +  g(aa+  +  CT-cf) 

+  B(t)dte~iut  +  B*{t)aeiuit,  (1) 

where  wq(t)  and  wr  are,  respectively,  the  qubit  and  the  res¬ 
onator  frequencies,  cr±  are  the  rasing  and  lowering  operators 
for  the  qubit,  a  (a  r)  is  the  annihilation  (creation)  operator  for 
the  resonator  photons,  g  (assumed  real)  is  the  qubit-resonator 
coupling,  B(t)  and  u>  are  the  effective  amplitude  and  the  fre¬ 
quency  of  the  microwave  drive,  respectively.  In  this  work  we 
assume  uj  =  ojr. 

For  the  microwave  drive  B(t)  and  the  qubit  frequency  ujq  (f ) 
[Fig.  1(b)]  we  use  Gaussian-smoothed  step-functions:  B(t)  = 
0.5i?o{Erf[(f  —  fB)/v/2o’B]  -  Erf [(f-fB  ~tb)/V2 ctb]}  and 
wq(f)  =  w0  +  0.5(Ao  -  A){Erf[(f  -  tq)/y/2crq]  -  Erf[(t  - 
fqeVv^Oqe]},  where  ts,  ts  +  Tb,  tq,  and  fqe  are  the  centers 
of  the  front/end  ramps,  and  <7b,  crq,  and  erqe  are  the  corre¬ 
sponding  standard  deviations.  In  numerical  simulations  we 
use  (7b  =  (7 qe  =  1  ns  (typical  experimental  value  for  a  short 
ramp)  while  we  use  longer  crq  to  make  the  qubit  front  ramp 
more  adiabatic.  Other  fixed  parameters  are:  g/2n  =  30  MHz, 
tb  =  1  ns,  ts  =  3  ns,  wr/27r  =  7  GHz,  and  wo/27r  =  6  GHz, 
so  that  initial/final  detuning  Ao  =  wr  —  ujq  is  1  GHz,  while  the 
disperse-stage  detuning  A  is  varied.  The  measurement  starts 
at  t  =  0  and  ends  at  if  =  fqe  +  2aqe  when  the  field  is  quickly 
released  [23]. 

Simplified  analysis- Let  us  first  consider  a  simple  dis¬ 
persive  scenario  at  large  qubit-resonator  detuning,  |A|  34> 
gy/n  +  1,  where  n  is  the  average  number  of  photons  in 
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FIG.  2:  (color  online),  (a)  Evolution  in  time  of  the  effective  field  am¬ 
plitude  Aeff  on  the  complex  phase  plane  for  initial  qubit  states  |0)  and 
1 1) ,  computed  numerically.  The  horizontal  line  emphasizes  slight 
asymmetry.  The  dots  indicate  time  moments  t  =  0,  5,  10,  15,  20, 
25,  and  30  ns.  (b)  Corresponding  probability  distributions  Pq(xv) 
and  P\{x<p)  for  measurement  (at  t  =  tf)  of  the  optimum  quadrature 
x,p.  Side  bumps  of  Po  and  Pi  are  due  to  non-adiabaticity.  We  used 
g/2n  =  30  MHz,  A/2n  =  50  MHz,  |Am|2  =9  [B0/2n  =  497.4 
MHz,  tb  —  1  ns),  (jq  =  3  ns,  fq  =  3.25  ns,  fqe  =  30  ns,  ts  =  3  ns, 
and  tf  =  32  ns. 


the  resonator.  In  this  case,  the  system  is  described  by  the 
usual  dispersive  Hamiltonian  [7]  Hd  =  (wo  +  g2/A)cr2/2  + 
(uv  +  crzg2  /  A)a^  a,  where  az  is  the  Pauli  matrix.  Af¬ 
ter  the  short  “catch”  stage  the  system  is  in  a  product  state 
(a|0)  +  /3|  1)) | Ain),  where  a  and  /3  are  the  initial  qubit  state 
amplitudes  and  Am  is  the  amplitude  of  the  coherent  resonator 
field,  Ain  =  —  i  J  B(t)dt  (so  n  =  |  Afn  | 2).  Then  during  the 
“disperse”  stage  the  qubit-resonator  state  becomes  entangled, 

ck|0)|Ao (t))  +  /3|l)|Ai(f)),  with  A0  =  Aine_^,Ai  =  Aine^, 

and  <j){t)  =  /g  [g2/ A(t')\dt'. 

The  distinguishablity  of  the  two  resonator  states  depends  on 
their  separation  |<5A|  =  |Ai  —  Ao|  =  2 1 Am |  sin  \<j>\  (see  numer¬ 
ical  results  in  Fig.  2).  The  released  coherent  states  are  mea¬ 
sured  via  the  homodyne  detection  using  the  optimal  quadra¬ 
ture  connecting  Ao  and  Ai,  i.e.  corresponding  to  the  angle 
(p  =  arg(Ai  —  Ao).  We  rescale  the  measurement  results  to 
the  dimensionless  field  quadrature  xv  =  {ae~lv>  +  a'elv>)/ 2, 
which  corresponds  to  the  </?-angle  axis  in  the  phase  space  of 
Fig.  2(a).  In  resolving  the  two  coherent  states,  we  are  essen¬ 
tially  distinguishing  two  Gaussian  probability  distributions, 
Paixtp)  and  P\{xv),  centered  at  ± |<5A|crcoh  with  crcoh  =  1/2 
being  the  coherent-state  width  (standard  deviation)  for  both 
distributions.  Then  the  measurement  error  has  a  simple  form 


i  r°° 

E=-  min(Po,  -Pi)  dx^ 

^  J  —  OO 


l-Erf(|gA|y^72) 


where  rj  =  t/coit/amp  is  the  detection  efficiency  [24],  which 
includes  the  collection  efficiency  r/t„i  and  quantum  efficiency 
of  the  amplifier  T/;llnp.  Unless  mentioned  otherwise,  we  as¬ 
sume  77  =  1,  which  corresponds  to  a  quantum-limited  phase- 
sensitive  amplifier  (for  a  phase-preserving  amplifier  77  <  1/2). 

Full  analysis.-ln  general  the  JC  qubit-resonator  interac¬ 
tion  (1)  is  non-linear  for  |Ain|2  >  A2/4 g2  [7]  and  the  res¬ 
onator  states  are  not  coherent.  The  measurement  error  E  is 
still  given  by  the  first  part  of  Eq.  (2),  while  the  probabil¬ 
ity  distributions  Po, i{xv)  of  the  measurement  result  for  the 
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qubit  starting  in  either  state  |0)  or  |1)  can  be  calculated  in 
the  following  way.  Assuming  an  instantaneous  release  of  the 
field,  we  are  essentially  measuring  the  operator  xv.  There¬ 
fore  the  probability  P{ xv)  for  the  ideal  detection  (rj  =  1) 
can  be  calculated  by  converting  the  Fock-space  density  ma¬ 
trix  pnrn  describing  the  resonator  field,  into  the  x  (,-basis,  thus 
obtaining  P{xv)  =  ^n{xv)pnm{t)'ilj^n(xlp)e-^n-rn'>(p, 
where  ij)n  ( x )  is  the  standard  nth-level  wave  function  of  a  har¬ 
monic  oscillator.  For  a  non-instantaneous  release  of  the  mi¬ 
crowave  field  the  calculation  of  P(xv)  is  non-trivial;  how¬ 
ever,  since  the  qubit  is  already  essentially  decoupled  from 
the  resonator,  the  above  result  for  P( xv)  remains  the  same 
[25]  for  optimal  time-weighting  of  the  signal.  In  the  case  of  a 
non-ideal  detection  (77  <  1)  we  should  take  a  convolution  of 
the  ideal  P{xv)  with  the  Gaussian  of  width  \J rj-1  —  1  crcoh. 
Calculation  of  the  optimum  phase  angle  p  minimizing  the  er¬ 
ror  is  non-trivial  in  the  general  case.  For  simplicity  we  still 
use  the  natural  choice  ip  =  arg(Aeff,i  —  Aeff,o),  where  the 
effective  amplitude  of  the  resonator  field  [26]  is  defined  by 
Aeff  =  n  V^Pn,n- 1-  The  field  density  matrix  pnm  is  calcu¬ 
lated  numerically  using  the  Hamiltonian  (1)  and  then  tracing 
over  the  qubit. 

Extensive  numerical  simulations  allowed  us  to  identify  two 
main  contributions  to  the  measurement  error  E  in  our  scheme. 
The  first  contribution  is  due  to  the  insufficient  separation  of 
the  final  resonator  states  |Aeff.i)  and  |Aeff  o}>  as  described 
above.  However,  there  are  two  important  differences  from 
the  simplified  analysis:  the  JC  nonlinearity  may  dramatically 
change  |<$A|  and  it  also  produces  a  self-developing  squeez¬ 
ing  of  the  resonator  states  in  the  quadrature  xv,  significantly 
decreasing  the  error  compared  with  Eq.  (2)  (both  effects  are 
discussed  in  more  detail  later).  The  second  contribution  to 
the  measurement  error  is  due  to  the  nonadiabaticity  of  the 
front  ramp  of  the  qubit  pulse  u>q(t),  which  leads  to  the  pop¬ 
ulation  of  “wrong”  levels  in  the  eigenbasis.  This  gives  rise 
to  the  side  peaks  (“bumps”)  in  the  probability  distributions 
Po,i(xv),  as  can  be  seen  in  Fig.  2(b)  (notice  their  similarity  to 
the  experimental  results  [9,  10],  though  the  mechanism  is  dif¬ 
ferent).  During  the  dispersion  stage  these  bumps  move  in  the 
“wrong”  direction,  halting  the  exponential  decrease  in  the  er¬ 
ror,  and  thus  causing  the  error  to  saturate.  The  nonadiabaticity 
at  the  rear  ramp  of  u>q{t)  is  not  important  because  the  moving 
bumps  do  not  have  enough  time  to  develop.  Therefore  the  rear 
ramp  can  be  steep,  while  the  front  ramp  should  be  sufficiently 
smooth  [Fig.  1(a)]  to  minimize  the  error. 

Now  let  us  discuss  the  effect  of  nonlinearity  (when  |  Ai„  | 2  > 
AW)  on  the  evolution  of  Aeff,o  and  Aeff,i  during  the 
disperse  stage.  Since  the  RF  drive  is  turned  off,  the  in¬ 
teraction  described  by  the  Hamiltonian  (1)  occurs  only  be¬ 
tween  the  pairs  of  states  |0,n)  and  |1  ,n  —  1)  of  the  JC 
ladder.  Therefore,  if  the  front  ramp  of  the  qubit  pulse 
is  adiabatic,  the  pairs  of  the  JC  eigenstates  evolve  only 
by  accumulating  their  respective  phases  while  maintain¬ 
ing  their  populations.  Then  for  the  qubit  initial  state  |0), 
the  qubit-resonator  wavefunction  evolves  approximately  as 
\ipo(t))  ~  e-Aml2/2  ^n(A?Ji/\/nT)e_*^o’n(t)|0 ,n),  where 


FIG.  3:  (color  online).  Optimized  measurement  error  E  vs  measure¬ 
ment  time  tf  (optimization  is  performed  over  A,  <rq,  and  iq).  The  mi¬ 
crowave  pulses  correspond  to  mean  photon  number  |Am|2  =6, 9, 12, 
and  15.  The  lines  are  guides  for  the  eye. 


the  overbar  denotes  the  (dressed)  eigenstate  and  (f>oin(t)  = 
Jt*  dt'[y/ A(i')2  +  4g2n— A(f')]/2  is  the  accumulated  phase, 
with  t d  =  £b  +  Tb/2  being  the  center  of  the  f?(£)-pulse, 
which  is  crudely  the  start  of  the  dispersion.  Similarly,  if 
the  qubit  starts  in  state  |1)  (following  the  ideology  of  Ref. 
[27],  we  then  use  1 10)  as  the  initial  state),  the  state  evolves 
as  | ipi(t))  ~  e-lAi”l2/2X:„(Ar„/V^!)e^i.’*(t)|l,n),  where 
</> i,n(t)  =  ftD  dt'[y/A(t')*  +  4 g2(n  +  1)  -  A(f')]/2.  Us¬ 
ing  the  above  definition  of  Aeff  and  assuming  |Ain|2  >  lwe 
derive  an  approximate  formula 


Aeff,o 


Ain  exp 


_  g  (]£ 

\J  A(F)2  +  4g2|Ain|2 


(3) 


The  corresponding  expression  for  Aeff,i  can  be  obtained  by  re¬ 
placing  —i  with  i  and  |  Ain|2  with  |  Ain  | 2  +  1.  These  formulas 
agree  well  with  our  numerical  results.  In  particular,  they  ex¬ 
plain  why  Aeff.o  rotates  slightly  faster  than  Aeff,i,  as  seen  in 
Fig.  2(a). 

Equation  (3)  shows  that  a  decrease  in  detuning  leads  to  an 
increase  in  the  rotation  speed  of  Aeff.  However,  in  the  strongly 
nonlinear  regime  |Ajn|2  A2/4g2,  the  angular  speed  sat¬ 
urates  at  d(arg( Aeff, o/i ))/dt  =  =Fg/2|Ajn|.  Thus,  the  rate  at 
which  the  Aeff,  i  and  Aeff,o  separate  is  limited  by 


d\6X\/dt  <  |g|,  (4) 

which  does  not  depend  on  |Ain|.  This  means  that  the  mea¬ 
surement  time  should  not  improve  much  with  increasing  the 
mean  number  of  photons  |Ajn|2  in  the  resonator,  as  long  as  it 
is  sufficient  for  distinguishing  the  states  with  a  desired  fidelity 
(crudely,  |Ai„|2  >  7/?/  for  E  <  10-4). 

Results  of  numerical  optimization.- Figure  3  shows  the  re¬ 
sults  of  a  three -parameter  optimization  of  the  measurement 
error  E  for  several  values  of  the  average  number  of  photons 
in  the  resonator,  |A;n|2  (assuming  77  =  1).  The  optimization 
parameters  are  the  qubit-resonator  detuning  A,  the  width  aq, 
and  the  center  tq  of  the  qubit  front  ramp.  We  see  that  for  9 
photons  in  the  resonator  the  error  of  10-4  can  be  achieved 
with  30  ns  measurement  duration,  excluding  time  to  release 
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FIG.  4:  (color  online).  Optimized  error  vs  measurement  time  if  for 
| Am | 2  =  9  and  quantum  efficiencies  rj  =  1  or  rj  =  1/2  (e.g.  for  a 
phase-preserving  amplifier),  taking  into  account  the  qubit  level  |2) 
(with  anharmonicity  A/2n  =  200  MHz)  or  assuming  a  two-level 
qubit. 


and  measure  the  field.  The  optimum  parameters  in  this  case 
are:  A/2ir  =  60  MHz,  aq  ==  4.20  ns,  and  tq  =  3.25  ns  [this 
is  a  strongly  nonlinear  regime:  |Ajn|2/(A2/4g2)  =9].  As 
expected  from  the  above  discussion,  increasing  the  mean  pho¬ 
ton  number  to  12  and  15  shortens  the  measurement  time  only 
slightly  (by  1  ns  and  2  ns,  keeping  the  same  error).  The  dotted 
blue  curve  in  Fig.  4  shows  the  optimized  error  for  |Ain|2  =  9 
and  imperfect  quantum  efficiency  p  =  1/2.  As  we  see,  the 
measurement  time  for  the  error  level  of  10-4  increases  to  40 
ns,  while  the  error  of  10~3  is  achieved  at  if  =  32  ns. 

So  far,  we  considered  the  two-level  model  for  the  qubit. 
However,  real  superconducting  qubits  are  only  slightly  anhar- 
monic  oscillators,  so  the  effect  of  the  next  excited  level  |2) 
is  often  important.  It  is  straightforward  to  include  the  level 
1 2)  into  the  Hamiltonian  (1)  by  replacing  its  first  term  with 
ccq 1 1) ( 1 1  +  (2wq  —  Al)|2)(2|,  where  A  is  the  anharmonicity. 
The  dispersion  can  then  be  understood  as  due  to  repulsion 
of  three  eigenstates:  |0,n),  |l,n  —  1),  and  |2,n  —  2).  As 
the  result,  Aeffp  rotates  on  the  phase  plane  faster  than  in  the 
two-level  approximation,  while  Aeff,i  rotates  slower  (some¬ 
times  even  in  the  opposite  direction).  The  Supplemental  Ma¬ 
terial  [28]  illustrates  evolution  of  the  resonator  Wigner  func¬ 
tion  corresponding  to  initial  qubit  states  |0)  and  |1).  The  de¬ 
tailed  description  of  the  effect  of  level  |2)  is  beyond  the  scope 
of  this  paper.  Here  we  only  show  the  optimized  error  for 
A/2tt  =  200  MHz  (a  typical  value  for  transmon  and  phase 
qubits),  |  Ain  | 2  =  9  and  rj  =  1  or  77  =  1/2  as  the  red  lines  in 
Fig.  4.  For  the  error  of  10“  3  the  measurement  time  if  is  31  ns 
and  39  ns,  respectively. 

The  presented  results  are  for  the  coupling  g/2n  =  30  MHz. 
While  the  dispersion  rate  scales  as  g  in  the  strongly  nonlinear 
regime  [see  Eq.  (4)],  the  scaling  of  the  overall  duration  of  the 
measurement  process  is  slower  than  g-1  because  of  signifi¬ 
cant  time  needed  at  the  initial  and  final  stages  [see  Fig.  1(b)]. 

In  this  work  we  are  not  interested  in  how  close  our  read¬ 
out  is  to  the  quantum  non-demolition  (QND)  measurement 
[29].  Notice  that  in  the  proposed  implementation  of  the  sur¬ 
face  code  [15]  the  measured  qubits  are  reset,  so  the  QND- 
ness  is  not  expected  to  be  important.  For  the  results  presented 


FIG.  5:  (color  online),  (a)  Time-evolution  of  the  quadrature  squeez¬ 
ing  (the  qubit  is  initially  in  state  |0}).  (b)  Measurement  error  vs  |<5A| 
obtained  numerically  (solid  lines)  and  using  Eq.  (2)  (dashed  line); 
the  evolution  stops  at  98  ns.  Here  g/2n  =  30  MHz,  |Ai„|2  =  9, 
aq  =  4ns,  and  tq  =  3.25  ns. 


in  Fig.  3  the  non-QNDness  (probability  that  the  initial  states 
|00)  and  1 10)  are  changed  after  the  procedure,  including  the 
release  stage)  is  crudely  about  5%,  which  is  mainly  due  to 
non-adiabaticity  of  the  rear  ramp  of  the  qubit  pulse.  It  is  pos¬ 
sible  to  decrease  the  non-QNDness  significantly  by  using  a 
larger  rear  ramp  width  crqe  and  correspondingly  increasing 
the  overall  duration  by  several  nanoseconds.  However,  the 
non-QNDness  cannot  be  reduced  below  the  level  of  few  times 
(g/A())2,  essentially  because  of  the  Purcell  effect  during  the 
release  stage. 

Squeezing.-The  JC  nonlinearity  causes  the  quadrature 
squeezing  of  the  microwave  field.  To  quantify  the  squeez¬ 
ing,  we  calculate  the  variance  Ax2  =  (a;2 }  —  ( x <p)2,  which 
can  be  put  in  the  form  As2  =  1/4  +  {a)a)/2  —  |(a}|2/2  + 
Re[((a2)  —  {a)2)e~2ltp]/ 2.  For  a  coherent  field  Ax2  =  1/4, 
thus  the  state  is  squeezed  [26]  when  4Ax2  <  1.  Notice  that 
the  degree  of  squeezing  depends  on  the  choice  of  tp.  However, 
in  order  to  see  the  effect  of  squeezing  on  measurement  error, 
we  compute  the  squeezing  along  the  measurement  direction, 
ip  =  arg(Aeff,i  —  Aeff,o)- 

Figure  5(a)  shows  evolution  of  the  squeezing  parameter 
4 Ax2  when  the  initial  qubit  state  is  |0)  (a  similar  result  is 
obtained  for  the  initial  state  1 1);  we  again  use  the  two-level 
approximation  and  r/  =  1).  Notice  that  at  first  the  field 
stays  coherent,  which  is  due  to  the  linearity  of  the  qubit- 
resonator  interaction  at  large  detuning.  Later  on,  however,  the 
interaction  becomes  nonlinear  due  to  decreased  detuning  and 
leads  to  quadrature  squeezing  reaching  the  level  of  ~50%  for 
A/27 r  <  100  MHz  (see  [28]  for  the  Wigner  function  evo¬ 
lution).  Figure  5(b)  shows  the  measurement  error  as  a  func- 
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tion  of  | <5A |  calculated  numerically  (solid  curves)  and  analyti¬ 
cally,  based  on  Eq.  (2)  (dashed  curve).  As  expected,  with  the 
squeezing  developing,  the  error  becomes  significantly  smaller 
than  the  analytical  prediction,  for  instance,  up  to  a  factor  of 
30  for  A/27T  =  250  MHz.  Notice  also  that  the  error  shown  in 
Fig.  5(b)  saturates  in  spite  of  increasing  separation  |<5A|.  This 
is  because  of  the  non-adiabatic  error  discussed  above. 

Conclusion.- We  analyzed  a  fast  high-fidelity  readout  for 
superconducting  qubits  in  a  cQED  architecture  using  the  con¬ 
trolled  catch,  dispersion,  and  release  of  the  microwave  pho¬ 
tons.  This  readout  uses  a  tunable  coupler  to  decouple  the  res¬ 
onator  from  the  transmission  line  during  the  dispersion  stage 
of  the  measurement,  thus  avoiding  the  Purcell  effect.  Our  ap¬ 
proach  may  also  be  used  as  a  new  tool  to  beat  the  standard 
quantum  limit  via  self-developing  field  squeezing,  directly 
measurable  using  the  state-of-the-art  parametric  amplifiers. 
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